Space positronium detection by radio measurements 
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Abstract 

The possibility of positronium detection at radio wavelengths is investigated in detail. 
For orthopositronium, the fine structure lines of the n = 2 level (1.62 cm, 2.30 cm and 3.48 
cm) may be probably detected by radiotelescopes of next generations with the sensitivity 
around \[iJy. The spin-flip line of the ground state (0.147 cm) may be observed from sources 
with a narrow annihilation line by modern telescopes with the sensitivity around lOfiJy. In 
the fine structure and spin-flip lines a maser effect is produced. The expected radio fluxes in 
the spin-flip line can reach detected limits for certain source geometries and optimal physical 
conditions inside the source. 
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Introduction 

The development of X-ray and 7-ray astronomy over the last ten years has shown that the 
(e + e~) annihilation process is exceptionally important in modern astrophysics; the presence 
in space of processes leading to creation of positrons has been unambiguously confirmed. 
However, a narrow 511 KeV line has only been detected in the Galactic center region (Tee- 
garden, 1994; Cheng et al., 1997; Purcell et al., 1997; Teegarden et al., 1997 and Smith, 
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Purcell, Leventhal, 1997). There are some sources, e.g. IE 1740.7-2942, Nova Muscae, and 
the Crab, whose spectra display bright transient emission features near 0.5 MeV (see the re- 
view of Smith et al., 1996). As theoretical analysis has shown, positron annihilation proceeds 
mainly via positronium states if the temperature of the annihilation medium is low (Gould, 
1989; Guessoum, Ramaty, Lingenfelter, 1991; Burdyuzha and Kauts, 1995 [hereafter BK]). 

The probability of positronium detection in the UV (A = 2431A) and IR (A = 1.31/xm, 3.75/xm) 
regions has been considered in detail in our previous articles (Burdyuzha, Kauts, and Yudin, 
1992; Burdyuzha, Kauts, and Wallyn, 1996). Here we consider the possibility of observing 
the radio line triplet at 8625 MHz (3.48 cm), 13010 MHz (2.30 cm), and 18496 MHz (1.62 
cm), corresponding to the fine structure of the orthopositronium n = 2 level and the spin-flip 
line at 203387 MHz (0.147 cm) of the hyperfine structure of the ground state. 

Since the observed narrow annihilation line sources are located in the Galactic Center 
region where the absorption at UV and IR wavelengths is high, the study of positronium 
detection at radio wavelengths is very important. 

Positronium atom 

The simplest positonium atom consists of an electron e~ and a positron e + , and can 
exist in two states - parapositronium with zero total spin of its electron and positron, and 
orthopositronium with the total spin equal to unity (see Fig. 1). The main distinction 
between a positronium and a hydrogen atom lies in the fact that in states with zero orbital 
momentum positronium is an unstable atomic system with lifetimes in the orthostate r ortho ~ 
1.33 x 10~ 7 n 3 sec and in the parastate r vara — 1.25 x 10~ 10 n 3 sec (here n is the principal 
quantum number). Also shown in Fig. 1 are the fine and hyperfine structure of the lower 
levels of positronium. As in the case of the hydrogen atom, the formation of positronium 
in an excited state leads to an electromagnetic cascade to lower levels and to appearence of 
characteristic lines (L a , Lp..., H a , Hp..., P a , Pp...). 

The lifetimes of lower states to annihilation are very short, so for low densities an equilib- 
rium population of levels has no time to form and populations are determined by processes 
of recombination or charge exchange followed by a radiative cascade and annihilation. In 
these conditions the appearence of maser radiation seems natural. 
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Positronium atoms are produced after the delay of positrons. The probability of positro- 
nium formation on one positron (/) depends essentially on the parameters of the medium (its 
temperature, the presence of dust, and so on). The value of / in different media has been es- 
timated by many authors (e.g. Bussard, Ramaty, and Drachman, 1979) and is ~ 0.9 4-0.95. 
For the case of an ionized hydrogen medium the rate of radiative capture in the process 
e + + e~ — > Ps + 7 is higher than the rate of direct annihilation e + + e~ — > 27 at temper- 
atures < 10 6 K (Gould, 1989). For the case of an atomic or molecular hydrogen medium 
an energetic positron, after entering the medium, effectively loses its energy by ionization 
and excitation of atoms and molecules of hydrogen and by Coulomb collisions. In this case 
annihilation processes proceed at low energies and positronium atoms will be formed at these 
energies by charge exchange reactions before annihilation. 

For astrophysical applications it is important that the formation of excited states of 
positronium be suppressed (BK). Therefore transitions between lower levels are more promis- 
ing for observations. 

Radio lines of the orthopositronium n=2 fine structure level 



The triplet of radio lines corresponding to n = 2 fine structure transitions may appear 
for radiative transitions 2 3 S\ — > 2 3 Pj (i = 0, 1,2) (see Fig 1). The lifetime of the 2 3 S\ level 
is determined by three-photon annihilation r(2 3 S'i) = 1.1/zs while the lifetimes of the 2 3 Pj 
levels are determined by the L a transition r(2 3 P) = 3.2ns. It is easy to see that a strong 
overpopulation of the 2 3 S' 1 level is created and that maser emission takes place in this triplet. 

The one-dimensional radiative transfer equation of in a line of this triplet is: 



dR hvi 



dx Auj 



A(2 3 S 1 -> 2 3 P)/4tt + B(2 3 S 1 -> 2 3 P i )Pjc 
A(2 3 S l Z 2 3 P t ) + E B(2 3 S 1 Z 2 3 P)P,Aiy C + 



annih. 



^ps{2 Pi) 



B{2 3 P l -> 2 3 S l )Pjc 
B(2 3 P t ^2 3 S 1 )P u AQ i /c + A Lo 



,i = 0,l,2 



where P v is the intensity in one of the lines of the triplet; n ps is the rate of positronium 
formation in the respective states per unit volume; A(2 3 S± — > 2 3 Pj) are the probabilities of 
spontaneous transitions 5.71 x 10 _6 s -1 , 5.96 x 10 -6 s _1 , and 2.89 x lO^s- 1 in levels 2 3 P , 2 3 Pi, 
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and 2 3 P 2 , respectively, and Bi are the Einstein coefficients. Afi; is the degree of background 
anisotropy; W^ nih is the probability of three-photon annihilation with the 2 3 Si level and 
A^ a is the probability of the L a transition. 

In this expression collisional processes are not included since densities in the annihilation 
region are low. For the case of ionized hydrogen this leads to the limit n<10 9 v / T c v cm -3 , 
where T is the Maxwell temperature of ionized hydrogen (Burdyuzha, Kauts, and Yudin, 
1992). For the case of atomic and molecular hydrogen our approximation will be valid over 
a wider range of densities. 

Since for any physical conditions the inequality A(2 3 Si — > 2 3 P)/47r <C B(2 3 S L — > 
2 3 Pj)4/c holds, we obtain: 



dx Aza- 



. B(2 3 S 1 -> 2 3 P l )Pjc 

1=0,1,2 



J B(2 3 P^2 3 ,S 1 )4A^/c + A La 



,i = 0,l,2 (2) 



For the conditions -B(2 3 Si — > 2 3 Pj)/*Af2j/c ~ ^Innih. a transition from the unsaturated 
to the saturated regime takes place. The respective brightness temperatures are T fe *' cr ' = 
(1.5, 1.0, andl.Z) x 10 11 (-^-) K. In the subsequent discussion we consider three possible 
cases. 

1) Tb < Tfr Cr (that is, the maser regime is not saturated). The probability of positronium 
production in the 2S, 2P states in different media has been considered many times (Nahar, 
1984; Hewitt et al, 1990; Biswas et al., 1991). Since A La > W^ nih and the formation 
rates h ps (2 3 Si) and n ps (2 3 Pi) differ insignificantly for positronium formation energies (BK), 
so *dff Sl) > " pj(23p<) and equation (2) takes the form: 



anmn. 



dP hvi 



dx Au. 
The corresponding depths Tj are: 



h ps (2 3 S 1 )—^B(2 3 S 1 - 2 3 P)/:/c 

' ' annih. 



hvt , - . 1 c 2 A(2 3 S 1 -> 2 3 P) , 

— -h ps (2 3 S 1 )— - —„ -l,i = 0,1,2 3 
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where / is a characteristic dimension of the annihilation region. 

Thus we see that r { ~ A^/uf ~ (2Jj + 1), that is r : T\ : r 2 = 1 : 3 : 5. 

As the simplest example, we consider a spherical annihilation line source radiating 
isotropically which has dimension I at a distance R from the observer. We take n ps (2 3 Si) ~ 
n ps (l 3 S'i)/8 (for the case of ionized hydrogen the asymptotic behaviour of these values was 
investigated by Burdyuzha et al. (1992); for atomic and molecular hydrogen the detail anal- 
ysis and the rationale for choosing this approximation was presented by BK). Then for the 
most intense line 2 3 S\ — > 2 3 P 2 iy 2 = 8625 MHz) the optical depth r 2 is 

T2 * 5 X 10_5 (^TS/) ( 10- Vi/cm 2 s) (lOlb) (li) (4) 



for Au/u ~ 1/200 (the width depending on the motion of positronium atoms (Teegarden, 
1994 and BK) and the intrinsic line width are of the same order). The value of ( g-iY/ ) ~ 
1 - 1.2 for / ~ 0.9 - 0.95. 

2) If Tft > T fe *' cr , the maser regime is saturated and its main characteristic will be a spectral 
flux density. In this case the ratios of the spectral flux densities depend on the values of 
B(2 3 Si — > 2 3 R i )F l/ Ail i . The total photon fluxes in the lines of the triplet are determined 
only by the rate of positronium formation in the 2 3 Si state: 

. = n(2 3 ^)/ 3 ^ BJlm { 4tt \ 

i 

where AQ F i is the anisotropy of the radio emission from the annihilation region in the 
respective line. As an example, for a flat spectrum and Af2j ~ constant we have: 

W 2»P 2 ) « 0.00! (j^JJf) {wS^i) (jfe) W 



3) For very high temperatures T b > T^ hi9h = (1.6, 3.3, and7.4) x 10 13 (^-) X determined 
by conditions B(2 3 Pi — > 2 3 Si)F v AQi/c m Al o , induced transitions dominate. Therefore the 
ratios of the spectral flux densities are F® : : ~ 1 : 3 : 5. For T fe > T^ h%9h line 
broadening due to induced transitions is the most important effect and an following increase 
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in temperature results in overlapping lines. The spectral flux densities in this regime will 
also be determined by formula (6) but it is necessary to take the increase in line widths into 
account adding the coefficient ~ (T b l,ht9h / T^) to this formula. 

In parapositronium the situation is more complex since the lifetimes of the 
levels r(2 1 Pi) ks t(2 1 S ) and accurate predictions require further study. 

The positronium spin-flip radio line 

The hyperfine structure transition l 3 ^ — > l 1 ^ of the ground state (spin-flip line) and 
the fine structure transitions can be observed in emission since the lifetime of the level 
1 3 <Si(t = 1.33 x 10~ 7 sec) is much longer than the lifetime of the l 1 ^ level (r = 1.25 x 
1CT 10 sec). Here the maser effect is also evident. This transition is of the Ml type with 
v = 203387 MHz, A = 3.37 x 1(T 8 s' 1 and intrinsic line width Av = 1.3 x 10 9 Hz. We 
have repeated the analysis for the spin-flip line similar to previous one. The transition from 
the unsaturated to the saturated regime occurs for B(l 3 Si — > l 1 S )I U AQ/ 'c ~ W^ nih (l 3 5i). 
The corresponding brigthness temperature is T b cr = 2.2 x 10 15 (^) K. The optical depth in 
the unsaturated regime for the conditions described in formula (4) is 

r ~ 5 x 10- ( ( ^ ( R ) '(J- V 2 f 7) 

\2 - 1.5 f J \10- 4 ph/ cm 2 s J \10 kpc J \l a.u. ) y! 

The spectral flux density in saturated regime from an annihilation region is 



The width depending on the motion of Ps atoms (Teegarden, 1994; BK) and the intrinsic 
line width are of the same order just as in the case of the fine structure lines. 

If T b > T^ i9h = 7.7 x 10 17 (|^) K, which is determined by the condition ^(l 1 ^ -> 
l^SijluAVl/c pa Wg^nfoX^So), then the spectral flux density in this regime will also be given 
by formula (8) but due to the increase of the line width the coefficient (3T£ i9h /4T b ) appears. 

Conclusions 
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Simultaneous observations of triplet radio lines can solve the problem of positronium 
identification. It follows from equation 4 that observations of such lines in the unsaturated 
regime are possible only from sufficiently compact and distant sources with the large flux 
in an annihilation line. Today such sources are unknown. From equation 6, the spectral 
flux densities in the saturated regime are just detectable with modern radio telescopes for 
( io-4ph/cm 2 s ) ( Aff 2 ) — Note that formula 6 gives the total spectral flux density from the 
entire annihilation region, and that the flux may be larger if the source is anisotropic. But 
it is certainly an unnatural proposal. 

For the spin-flip line the optical depth in the unsaturated regime is very small since 
this is the Ml type transition and it is impossible to detect. But the situation is different 
in the saturated regime. The spectral flux density in the saturated regime is higher than 
that in the triplet lines (equation 8) but the point of transition from the unsaturated to 
the saturated regime (T fe cr ) is shifted. As was shown by Kardashev (1979), this transition 
is optimal for interstellar communication but so far attempts to observe possible planetary 
systems at 203387 MHz have not been successful (Mauersberger et al., 1995). It is very 
desirable to extend these observations to possible sources of the 511 keV annihilation line 
with high brightness temperatures in the mm range (radio jets sources, radiopulsars and so 
on). The search for this radio line can be done in the direction of the Galactic center but 
to date a narrow annihilation line point source has not been detected there, so the exact 
predictions cannot be done. One promising source for spin-flip detection is radiopulsars. 
The problem is a difficult one, since necessary conditions are the presence of positrons (e.g., 
in a pulsar wind), and an appropriate deceleration length; a medium with a define density 
and dimension is required. 

In addition to radio observations positronium can be also detected in the recombination 
lines L a (243lA) , if Q (1.31/zm) and P q (3.75/xto). This possibility was investigated in our 
previous articles (Burdyuzha, Kauts, and Yudin, 1992; BK; Burdyuzha, Kauts, and Wallyn, 
1996). The L a line is difficult to observe in the direction of the Galactic center because of 
high absorption but it may be observed from active galactic nuclei with UV bumps in their 
spectra and from new point sources in our Galaxy which have to be detected. Observations 
of L q (2431t4) can be carried out with the HST. The IR lines H a (1.31fj,m), ifg(0.97//m) and 
P a (3.75fxm) can be observed by Keck telescope and probably by UKIRT, IRTF IR telescopes. 
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